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A particularly simple model belonging to a wide class of coupled maps which obey a 
local conservation law is studied. The phase structure of the system and the types of the 
phase transitions are determined. It is argued that the structure of the phase diagram 
is robust with respect to mild violations of the conservation law. Critical exponents 
possibly determining a new universality class are calculated for a set of independent 
order parameters. Numerical evidence is produced suggesting that the singularity in the 
density of Lyapunov exponents at A = is a reflection of the singularity in the density 
of Fourier modes (a "Van Hove" singularity) and disappears if the conservation law is 
^v' \ broken. Applicability of the Lyapunov dimension to the description of spatiotemporal 

chaos in a system with a conservation law is discussed. 
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I. INTRODUCTION 



Coupled map lattices (CML) with an additive conserved quantity became a subject of intensive 
^^ ' research recently pj-pl. On the one hand such CML's are often obtained as phenomenological models 

f^ . representing the dynamics of a large number of interacting macroscopic structures. On the other hand 

04 ' they are a natural result of finite-difference approximations of continuous nonlinear partial differential 

^^ , equations such as the Kuramoto-Sivashinky equation y, or a phenomenological Cahn-Hilliard equation 

B describing the nonlinear dynamics of several systems with conscrvcd-ordcr-paramctcr. 
Models of this class are expected to represent several typical non-equilibrium physical phenomena. 
^ • For instance, surface waves Q, where the average depth of the fluid in the container is conserved, elec- 

^p»' trohydrodynamic instabilities in nematics with insulating plates 0|, where the total charge is conserved, 

disturbances in the atmosphere and ocean systems, where the total (depth integrated) heat is conserved, 
Q ' and even some types of hard turbulence [p[ . As such they are significant as tools for studying the coni- 

es , plex spatiotemporal behavior of spatially extended nonlinear systems, especially in the strongly chaotic 

1-^ ■ regime, where the analytical methods designed for weak nonlinearities become ineffective. 

For coupled map systems with an additive conserved quantity several major points are still awaiting 
clarification. First of all, what is the effect of the conservation law on the structure of the phase 
diagram and the character of the phase transitions? What is the connection between the conservation 
law and singularities observed in the density of Lyapunov exponents p],p||lO| and what is the origin 
C^ ' and significance of these singularities? Another important issue is to determine which parameters best 

describe the dynamics of an extended system and what their limitations are. In particular, it is unclear 
whether the Lyapunov spectrum provides any exclusive information about the chaotic dynamics that 
cannot be obtained by other methods. And finally, we would like to know whether the reduction of the 
system dynamics to symbolic form preserves the main characteristics of the chaotic dynamics and can 
provide us with the complete description of the latter. 

The model chosen should be relatively simple yet represent most of the typical features under con- 
sideration. Most important, it should have nontrivial phase diagram. With this in mind we pick the 
1-dimensional collection of L diffusively coupled chaotic maps: 

<+' - < + (/«-i) ~ 2/«) + /«+i)) (1) 

with periodic boundary conditions imposed. The local map was chosen to be 

f{x) ~ ax + bz{l — z), z — frac{x). (2) 

This CML can also be regarded as a finite difference approximation of the differential equation con- 
tinuous in both space and time 



dtu{x,t)=dlf{u{x,t)). (3) 

A differential equation of this form represents the competition between two opposing tendencies: gen- 
eration of chaotic perturbations by the nonhnear part of f{x) and dissipation of these perturbations by 
the diffusive couphng introduced by the second spatial derivative. 
The model clearly possesses the conservation law: 

1 ^ 
u — Y z_^ ^^ ~ const, (4) 

so that aside from the two parameters of the local map a and b we have an additional control parameter 
- the additive conserved quantity u, which is defined by the initial condition and is of critical importance 
to the behavior of the system. 

We are primarily interested in the dynamics in the "thermodynamic limit" L ^ oo, although the 
numerics is obviously restricted to finite systems. 

The outline of this pap er is as follows. In section O we study the phase diagram of the coupled 



map model. In section [II we introduce the symbolic (reduced) description of the system dynamics. In 
section |^ we discuss the quantitative description of the reduced dynamics. In section M we study phase 
transitions in our model and determine their types. In particular, we study the effect of the conservation 
law on the type of the transition and on the values of critical exponents. In section N^ we discuss the 
applicability of the Lyapunov dimension to the description of the dynamics in the system. In section 



VII we present numerical data sugg estin g the reason for the existence of the singularity in the spectrum 



of Lyapunov exponents. In section VIII we demonstrate the effect of violations of the conservation law 



on the system dynamics. The paper ends with a summary and discussion in section [X 



II. PHASE DIAGRAM 

Despite its simple form, this model has a very rich structure. Numerical simulations show that depend- 
ing on the values of the control parameters it can be strongly or mildly chaotic, show spatiotemporal 
intermittency (STI), give rise to pattern formation or simply decay into the spatially uniform stable 
state, or (for L - even) a 2-cycle in both space and time. Both asymptotically regular (non-chaotic) 
states can be described by a single equation: 

u] ^u+{-iy+'"A. (5) 

In order to gain some insight into the phase diagram we analytically determine the boundaries of the 
stability regions of the two non-chaotic asymptotic states of the system. Linear stability analysis of the 
spatially uniform state gives Lyapunov exponents (equivalent to the growth rates) in the form: 

A„ = ln|l-4(a + 6-26u)sin2(^)|. (6) 

From this one can conclude that the region of linear stability of the uniform phase is given by 

a + b 1 a + b ,„, 

Note that at the upper boundary all exponents corresponding to all Fourier modes k = ^j^ become real 
and positive, whereas at the lower boundary the k = n mode goes unstable to a period 2 oscillation. 

Stability analysis of the square of the map (|^) gives the boundaries of the linear stability region of 
the 2-cycle state written in a slightly different form: 

u+ if i + n is even, ,„. 

u_ if z + n is odd. 

The dynamics is invariant under the transformation u — > ±1. As a consequence it turns out that 
there are two types of 2-cycles possible. One (type-I) with 

[u-] = [u+] - 1 (9) 



([•] denotes the integer part) requires 

26 , 

u±=u(l±- — (10) 

and has a stabihty region bounded by the surfaces given by the foUowing two equations: 

(1 - 2a + 26 + 46m_)(1 - 2a - 26 + 46m+) = -1 (11) 



and 



26- (1^) 



At the upper boundary dl3) the k — tt mode becomes growing, while at the lower boundary (11) a Hopf 
bifurcation of the |A:| = tt/2 modes occurs. This is quite fortunate, since one can obtain the analytic 
expressions for the phase boundaries for a system of arbitrary size L from the analysis of a system with 
L = 4. 

The other (type-II) 2-cycle is such that 



= u and [u_] = [w+] (13) 

and, for the local map f{x) given by eq. (0), it can only exist at the stability boundary of the uniform 
state given by 

a + b 1 

u = , (14) 

26 46' ^ ' 

but can have an arbitrary amplitude A — (w+ — w_)/2, subject only to the condition (l3h. 

Fig. |l| presents two cross sections of the parameter space. We will denote the region where the 
uniform state (1-cycle) is linearly stable as the phase LI. Similarly, the phase L2 will stand for the 
linear stability region of the 2-cycle state. As we are going to see later, the 2-cycle state is not the only 
possible asymptotic state in this phase, so it will be useful to introduce the additional subdivision of 
this phase into sub-phases for a more detailed analysis. 

The attractors of the phases Tl and T2 are chaotic. The two phases are not essentially different. 
One can easily find a continuous trajectory in the parameter space that would join arbitrary points in 
phase Tl with those in phase T2 without intersecting any phase boundary. Although the dynamics of 
the system is somewhat different in the two phases for the set of parameters used, this distinction is 
introduced mostly for convenience. 

An important property of the system is that despite the large number of degrees of freedom the type of 
the attractor (and therefore, the type of behavior, if one excludes long transients) seems to be uniquely 
determined by the values of control parameters a, 6 and u and is independent of the details of the 
initial state. If there exists a single attractor, then the basin of attraction is (almost) all configuration 
space. In this case averages over the attractor can be estimated as the time average from a single initial 
condition, and we can call the system ergodicQ. 

Indeed, numerical data suggest that the attractor is unique in most of the parameter space. However 
it is not always the case: close to the boundary Ll-Tl a frozen pattern may form and the details of a 
pattern do depend on the initial conditions. So, on the timescale used in our calculations (typically of 
order 1 million iterations) the system did not appear ergodic. Of course, that does not mean that the 
system could not become ergodic on a yet larger timescale. 

Another exception is the regions inside the ordered phases LI and L2, where two attractors, chaotic 
and non-chaotic, can coexist. We will later discuss this situation in more detail. 

One of the objectives of this paper is to study the effect of the conservation law on the dynamics of the 
extended chaotic system. We therefore would like to follow the changes in the parameters characterizing 



^We assume that, at least with respect to the numerical computation, there exists a "physical measure" such 
as discussed by Eckmann and Ruelle m_U that eliminates any ambiguity in the choice of invariant measures on 
the attractor. 



the dynamics as a function of the conserved quantity u. In further work we will fix the values of the 
other two parameters at a = 0.4, 6 = 1.3, which (as seen from the phase diagram, fig. |l|) will allow us 
to study the regimes of interest (various chaotic as well as periodic states). 

Looking at the phase diagram, one can expect that this model should experience at least four bifur- 
cations or phase transitionsn as u is varied in the interval < u < 1. Equations (O) and (O) give us 
Ua ~ 0.0520 and Ub — Ud~ 0.5 « 0.1538 as the boundaries of the 2-cycle stability region. According to 
eq. (0) the uniform state loses its stability at Uc ~ 0.4615 and Ud ~ 0.6538. 

Phase transitions from ordered to chaotic states are common occurrences in coupled map lattices [|l2| . 
They may be continuous or discontinuous. For continuous transitions we might expect the transitions 
to fall into various universality classes, characterized by scaling exponents for various diagnostics of the 
chaos near the transition. In this case the transitions in CML's may be representative of transitions 
to spatiotemporal chaos in more general extended non-equilibrium systems. Although the qualitative 
features determining the universality classes are not understood one expects that symmetries, such as 
the Ising symmetry studied by Miller and Huse ||l^, and conservation laws, rather than the detailed 
properties of the local maps, will be important. 

In particular, it has been suggested [l3 that under certain very general conditions (e.g. in systems 
with a unique absorbing state) the transition should fall into the universality class of directed percolation, 
although some counter examples to this statement are known. It is nevertheless interesting to check 
whether any of the phase transitions in our model belong to the universality class of directed percolation, 
since, as we are going to see below, the absorbing state in our model is in fact unique for any choice of 
control parameters. 



III. REDUCED DYNAMICS 

In order to understand the spatial dynamics of the system better and to see the finer details of the 
phase diagram we will (following Kaneko ||l5[) reduce the description of the dynamics to a finite number 
of states: in terms of this reduced dynamics each site of the lattice can be marked either "laminar" or 
"turbulent" thus making up a set of laminar and turbulent domains. Then one would naively expect 
a laminar domain to be a region of the lattice with a relatively slow chaotic dynamics (no large, if 
any, positive local Lyapunov exponents ||l6|] ); and a turbulent domain to be a region where the chaotic 
dynamics is fast (with at least a few large positive local Lyapunov exponents). Positive Lyapunov 
exponents will inevitably make the turbulent domains spatially irregular, while laminar domains tend 
to be spatially regular. 

We will not put the terms "laminar" and "turbulent" in quotes below, nevertheless one should clearly 
understand that these are just a convenient notation and thus are restricted in meaning. 

We need a simple criterion that will determine whether a given site belongs to a turbulent or laminar 
domain (we will only consider as laminar states those that are close to either uniform or 2-cycle con- 
figurations). The simple way to distinguish between (uniform) laminar and turbulent sites numerically 
would be to call a site j laminar on time step n if 



u 



ji < e and |u" - u"_^i| < e (15) 

and turbulent otherwise. The problem with this definition is that any 2-cycle with amplitude 

would be considered turbulent, which is clearly incorrect. Therefore, the definition of a laminar domain 
has to be generalized to include a zig-zag pattern with slowly varying envelope of arbitrary amplitude: 






< 6. (17) 



^ Since we are interested in diverging correlation lengths near the transitions between different states in the 
L — > oo limit, and the consequent possibility of universal exponents, we will use the term "phase transitions" 
rather than "bifurcations". 



In this particular model we set e ~ 0.01. 

Based on this reduction one can distinguish between various types of dynamics in the system. Pictures 
representing the spatiotemporal evolution of the system in symbolic form are so characteristic that one 
can easily determine which part of phase diagram the system is in just by looking at the patterns. 
We will examine several typical pictures, highlighting the most interesting phenomena observed in this 
model. 

The behavior of our CML in the phase L2 is nontrivial. We already know that inside the phase there 
is a stable ordered state. Our numerical data imply however that this ordered 2-cycle state is stable in 
the nonlinear sense only for values of u satisfying Ua < Up < u < Un < Ub (figure |2|). As we are going to 
see later, it is quite hard to determine the critical values Up and u„ exactly, but they seem to approach 
the fixed values Up ~ 0.063 and m„ « 0.082 in the thermodynamic limit L ~> cx). 

For Up < u < Un the ordered state is reached through a chaotic transient whose lifetime is usually 
very small, no matter what initial condition is chosen, and thus the asymptotic behavior dominates. For 
relatively small systems (L ^ 400) this results in a fast decay of any initial state into a limit 2-cycle. For 
larger systems the quiescent asymptotic state may be completely regular as well, but it may also have 
a few localized turbulent defects moving with unit speed through the homogeneous 2-cycle background. 
Defects moving in opposite directions eventually die out colliding with each other. Similar properties 
were observed in a number of 1- and 2-dimensional models featuring spatiotemporally intermittent 
dynamics (see |17| and references therein for example). 

For Ua < u < Up as well as for Un < u < Ub the ordered state is only conditionally stable (stable to 
small perturbations) and most initial conditions result in a spatiotemporally chaotic asymptotic state 
consisting of a combination of laminar and turbulent domains, with the laminar domains featuring 
exactly the same structure as the ordered state: the type-I 2-cycle. Therefore the phase L2 can be sub- 
divided into three sub-phases according to whether any turbulent domains are present in the asymptotic 
state together with the laminar background whose structure is the same throughout the phase L2. 

Fig. y shows schematically the (time averaged) fraction of the lattice occupied by turbulent domains 
in the phase L2 as a function of parameter u. In the sub-phase L2i (up < u < w„) the asymptotic state is 
laminar. In the sub-phases L2p {ua < u < Up) and i2„ (u„ < u < Ub) the asymptotic state can be either 
laminar or spatiotemporally chaotic. It is interesting to note that the characteristic patterns produced 
by the turbulent domains are different in the two sub-phases featuring a persistent spatiotemporally 
intermittent state. As a result one might expect to see two quite different phase transitions at u — Up 
and u = Un- 

A pattern typical for the sub-phase L2p is presented on fig. p|(a). This state is very similar to some of 
the spatiotemporally intermittent states of a (1-dimensional) model studied by Chate and Manneville 
p9| . The major difference is that in our model turbulent domains typically have a larger size. Since the 
absorbing (laminar) state is an ordered one and is unique, we might expect a phase transition at Up to 
belong to the universality class of directed percolation jlj] . 

In fact, this kind of STI state is not specific to discrete extended systems. A very similar state can 
be observed in some continuous models as well, for instance po| , in a damped Kuramoto-Sivashinsky 
equation: 

dtu{x, t) = —rju — d^u ~ d^u — udxU (18) 

with 77 w 0.075. 

Fig. H(b) presents another type of spatiotemporally intermittent state that can be observed in the 
phase L2. Its distinguishing feature is that it is composed of a set of virtually immobile turbulent nuclei 
"mediated" by creation and absorption of nonlinear waves, propagating trough the laminar background 
with unit velocity. This type of "nuclear" STI state is characteristic for the sub-phase L2„. 

The phase LI also features internal structure (see fig. ||). For Uc < u < uj (sub-phase Lli) any 
initial configuration decays quickly into a uniform stable state. In other words, the uniform state is 
nonlinearly stable. For Uf < u < Ud (sub-phase Llf) the uniform state is only conditionally stable 
and large deviations from it result in a spatiotemporally chaotic asymptotic state. The transition point 
separating the two sub-phases is estimated to be Uf « 0.53 in the thermodynamic limit. 

Fig. ||(d) represents a frozen pattern, characteristic of the chaotic asymptotic state of the system in 
the sub-phase Llf and in the phase Tl close to the boundary with LI (we will denote this region Tli). 
It is probably more appropriate to call this type of dynamics locked chaos: the chaotic state is almost 
stationary in terms of reduced dynamics, with chaos localized in turbulent domains. As was briefly 
mentioned in section O, the dimensions and locations of the turbulent domains depend on the details 
of the initial state of the system. The laminar state in this sub-phase is uniform. 



And finally, fig. ^(c) gives the reduced dynamics representation of critical behavior displayed by the 
CML in the chaotic phase T2 close to the boundary with LI. Here we encounter yet another example of 
a spatiotemporally intermittent state observed in our model. The laminar state is again a 2-cycle (more 
specifically the type-II 2-cycle) and it is also absorbing, i.e a new turbulent domain (we will call these 
defects because of their small size) can never originate inside a laminar domain, it can only be spawned 
by other defect (s). However a defect can be consumed by a laminar domain, or alternatively it can be 
destroyed in a collision with another defect. 

The most prominent feature of this picture is "spontaneous" creation and annihilation of turbulent 
pulses (defects) moving in different directions with different (but constant) velocities. Therefore we 
may alternatively regard these defects as traveling waves. Naively one would expect that the condition 
^ :^ 1, where ^ is the correlation length, is necessary and might also be sufficient for the formation of a 
number of traveling waves. Numerical results for our model support this assumption (in disagreement 
with the stronger restriction [pl| , according to which the correlation length ^ should be comparable to 
the size of the system L). 

Nevertheless, since in the strongly chaotic regime the correlation length is of order one lattice spacing, 
the condition ^ ^ 1 is usually only satisfied close to the hypersurfaces in the parameter space on which 
the correlation length diverges, i.e. where a continuous phase transition occurs. This is clearly the case 
of fig. H(c): we have a continuous phase transition at u = Uc ~ 0.4615. 

Deeper in the chaotic phases Tl and T2, away from the phase boundaries, strongly chaotic behavior 
could be observed. Here almost all sites on the lattice exhibit turbulent behavior, and only occasionally 
a laminar domain of a very small size is created and then quickly consumed by the neighboring turbulent 
sites. We would call this type of dynamics strong chaos in contrast to the mild chaos, where all chaotic 
dynamics is localized to turbulent domains, occupying only a part of the lattice, while the rest of it is 
in the laminar state. 



IV. DOMAIN LENGTHS 

A quantitative description of the reduced dynamics is provided by the probability distribution func- 
tions Pt{l) and Pi{l) giving the probability for a turbulent (laminar) domain to have length /. In order 
to calculate these functions numerically we used a single random initial condition and let the system 
evolve, counting how many times a laminar (turbulent) domain of a given size formed. The resulting 
distributions did not depend on a particular choice of the initial condition for any given set of control 
parameters corresponding to ergodic dynamics, i.e everywhere except the sub-phase Llf. 

Our calculations show that in the chaotic state (away from phase boundaries) both Pt{l) and Pi{l) 
decay exponentially. Even more important, typical lengthscales appear to be almost independent of the 
value of the conserved quantity u, being h fv 1 for laminar and It ~ 10^ for turbulent regions. This is 
consistent with the point of view that most of the system is in a turbulent state and the probability of 
encountering a laminar region at any particular location is very small (but finite and independent of the 
system size) and decreases rapidly with increasing domain length, while there is no spatial structures 
defining alternative lengthscales. 

The behavior of length distribution functions might be expected to change substantially as the system 
gets close to the phase transition points. For example, if the transition is continuous the correlation 
length may grow and the critical effects might introduce alternative lengthscales that would modify the 
form of the distribution functions. 

We first examine the "percolating" STI state characteristic of the sub-phase L2p. Fig. 0(a) shows 
that both Pt{l) and Pi{l) decay exponentially for large I. This is compatible with the assumption that 
the phase transition a,t u = Up might in fact belong to the universality class of directed percolation p8| . 
It is interesting to note though, that Pi has two branches, one corresponding to even, the other to odd 
size of a domain. About the only useful information that one can extract from this data is typical length 
scales of laminar and turbulent domains, which, for u = 0.06, appear to be tens of lattice spacings for 
both. 

In the previous section we saw that the STI states observed inside L2p and i2„ are considerably dif- 
ferent. As a result fig. §(b), which corresponds to the "nuclear" state, differs from fig. ^a) substantially: 
though Pt{l) still decays exponentially, Pi{l) does not, but has another peak at Is « 200. This peak 
is not a finite size effect (which can be shown using a larger system) and indicates the presence of an 
internal spatial structure with characteristic length Is (typical separation between the "nuclei") in the 



STI state. The typical width of the "nuclei" is, in turn, determined by Pt{l)- 

Numerical data for a larger system (L = 4096) suggest (in contrast with the results obtained by Chate 
and Manneville [^ ) a crossover type of behavior for the distribution of laminar domain lengths away 
from the onset of STI: 

,. JZ", a<0 for !<?<?„ , 

^^^^'^\exp{-l/Q for/>/,. ^^^> 

As a result, if the characteristic length Zs — » oo as w — > m+ we should expect a pure power law decay of 

Pl{l) at U = Un- 

Another interesting phenomenon, can be pointed out in fig. ^(c). Both Pi{l) and Pt{l) have two 
branches, one corresponding to even, the other to odd length of a domain. The behavior of these 
branches is quite peculiar, they cross at some crossover length Icr ~ 36: 

peven^j^ > P^'^'\l), I < l^r] (20) 

Pf'"="(/) < P{"^'^{1), I > Icr- (21) 

The difference between the branches is larger for L-even and smaller for i-odd. Similar relations hold 
for Pt{l), though the difference between the branches is less pronounced than in the case of Pi. 

Close to the boundary T2-L1 a typical state (fig. H(c)) is composed of a collection of laminar domains 
separated by small turbulent defects. As we are going to see later, most of the turbulent domains tend 
to have a fixed length It = 6, independent of the distance (in parameter space) to the transition point. 

All defects move with a constant velocity, but while the majority of defects is moving with the maximal 
speed V = ±1 the rest have a smaller speed \v\ < 1. Therefore we have a considerable probability of 
encountering a laminar domain bounded by a pair of defects which move with equal and opposite speeds. 
If at some particular time n the length of such a laminar domain was even (odd), it will remain even 
(odd) for as long as these defects exist since the length will increase (decrease) by 2 at each time step. 

It is reasonable to assume that the details of the "defect interaction" favor the creation of domains 
with length of a given parity, say even over odd. As a result, one will see that the probability of finding 
a laminar domain of small and even size I is higher than the probability of finding a laminar domain 
with comparable odd size / — 1 (or I + 1), thus splitting the function Pi{l) into two branches. 

Since chaotic fluctuations tend to destroy the deterministic predictions like the one we just discussed, 
we cannot make any rigorous conclusions about the dynamics of larger laminar domains. The numerical 
data suggests that the characteristic lengthscales determining the decay rates of the two branches are 
different, so they cross at some crossover length Icr determined by the average distances between the 
defects moving with different speeds. 

In the case of a frozen pattern one should not expect the probability density to decay exponentially, 
rather it should display several peaks, broadened due to the chaotic fluctuations of the domain bound- 
aries, at the lengths occurring in a particular pattern and then go sharply to zero. This applies to both 
Pt and Pi. Fig. 0(d) satisfies this prediction quite well. The underlying reason of this kind of behavior 
is the non-ergodicity of the system in the frozen pattern forming regime. 

The numerical results in fig. imply that in the model studied in this paper simply the ergodicity 
of the system dynamics is a sufficient condition for the probability distribution function Pt{l) to decay 
exponentially at lengths greater than the correlation length f in any disordered state. On the contrary, 
the distribution of laminar domains is more specific and informative. Pi{l) still decays exponentially in 
the strongly chaotic states. In the spatiotcmporally intermittent state however the behavior of Pi varies 
widely: it might or might not decay exponentially. For instance, in case of the "nuclear" STI state we 
observe the crossover from exponential to power law type of decay. Similarly to |19| we expect to see a 
pure power law decay at the transition point u = u„. 



V. PHASE TRANSITIONS 

A. Order Parameters 

As mentioned above, we expect to have 4 distinct phase transitions, in which the system goes from 
either uniform state or 2-cycle to a chaotic state. The transition to chaos in this example of an extended 
system with a local conservation law does not follow the period doubling cascade or other routes to 



chaos characteristic of low dimensional dynamical systems, instead the system goes directly from a 
simple dynamical state (fixed point, period-2) to a chaotic state. In other words it has a character 
similar to phase transitions in Hamiltonian statistical systems where a symmetry of a basic state is 
destroyed upon crossing of the critical point. This feature is common to all the phase transitions in this 
model. 

A conventional dynamical systems approach to the treatment of phase transitions in a deterministic 
chaotic system would be to calculate the maximal Lyapunov exponent, Xmax- The bifurcation from the 
ordered to disordered state then occurs at the values of the parameters where the exponent changes 
sign and becomes positive. In the case of CML with a conservation law, one of the exponents is always 
zero, therefore Xmax does not change sign, but increases (continuously or discontinuously, depending on 
the type of transition) from zero as the system crosses the boundary between ordered and disordered 
phase. The maximal Lyapunov exponent can be considered as an example of a global (intensive) order 
parameter. 

On the other hand, in order to get some additional insight into the spatial dynamics of the system, 
it might be advantageous to introduce a local (extensive) order parameter. A good candidate seems to 
be the density h of the Kolmogorov-Sinai entropy, Sks- 

h=yKS = 2 ^ ^"- ^^^^ 

Am>0 

It is clearly zero in the ordered phase and positive in the disordered phase. 

The calculation of the Lyapunov spectra is numerically costly, so a different approach is used for 
comparison, based on the reduced description of the system's dynamics. In terms of reduced dynamics 
the phase transition from the ordered to disordered state can be represented as the appearance of a 
set of disjoint turbulent domains on the laminar background. This leads us naturally to the measure 
of the set of turbulent domains, pt, as another local order parameter, describing how the laminar state 
becomes turbulent. In the ordered state pt — 0, in the disordered state < pt < 1. 

And, finally, let us introduce yet another order parameter, Cch- One can decompose the chaotic 
dynamics into modes using the Karhunen-Loeve decomposition p2[ . The mode intensity Ek is defined 
as the eigenvalue of the integral equation 

L 

Y,K{j,i)Mi) = EkMJ), (23) 

1=1 

where the kernel 

i^(z,j) =<<"">« (24) 

is just the 2-point correlation function. Due to the translational invariance of the system (in case of 
ergodic dynamics) K(i,j) = C(i — j) and therefore the eigenfunctions V'fe(*) ^^'6 just Fourier modes. 
Consequently, the eigenvalues Ek are given by the values of the static structure function 

Ek = S{k) =< |4f >„, (25) 

where u^ is the Fourier transform of the map variable. 

The total intensity of the dynamics is defined as E = ^^ Em ■ It includes the contributions from the 
chaotic as well the non-chaotic modes. 

In the non-chaotic phases LI and L2 the stable stationary state is given by the general formula (ph. 
Therefore, the structure function might only be nonzero at two values of the wave vector k: 

5(0) = M^ (26) 

and, for a zig-zag state, 

S{tt) = A^. (27) 

In other words, if we want the order parameter to represent the strength of chaos in the system, 
it should be defined through the intensity of the chaotic modes only. So we arrive at the following 
expression 



T-1 

1 v-^ , , 27rTO . 

ech = J ^ o{k,n), kjn = —j—, (28) 

m—l 

which is identically equal to zero in any ordered phase and larger than zero in any disordered phase. 

All order parameters introduced are expected to become asymptotically independent of the system 
size L in the thermodynamic limit L — > cx). 



B. Continuous transition 

Most of the attention in this study was devoted to the order-disorder transition occurring at the 
boundary L1-T2. This transition is easy to study and it is expected to be quite common in models 
described by a CML with diffusive coupling and conserved map variable density. 

The transition point is defined by the equation f'{uc) = 1/2. One can see from eq. (|g) that a.t u = Uc 
a period 2 bifurcation occurs and the k = tt mode becomes growing, making the uniform configuration 
unstable. 

Details of the transition (see next section) suggest that this transition is a continuous (second order) 
phase transition, and the numerical data support this conclusion. One of the clear indications of this 
fact is presented in fig. o, which shows that the correlation length diverges as we approach the phase 
transition point 

^ ex (uc - u)-" (29) 

with critical exponent v estimated to be of order 0.8. The correlation length is hard to measure however, 
and the precision of this result is low, so that this value cannot be considered reliable. 

A better diagnostic of a diverging lengthscale is the average length of a laminar domain k = J^i lPi{l) 
which is seen to scale (fig. ^(a)) as 

li(x{uc~u)-^', /i« 1.0 ±0.02 (30) 

which implies, in particular, that very close to the transition point the lattice configuration consists of 
a few large laminar domains, separated by turbulent defects of finite size. 

The average length of a turbulent domain It — ^ilPtil) does not scale (fig. ^(b)), but shows an 
exponential dependence on the distance from the critical point: 

k ^ IT exp{— ), ut » 0.006, (31) 

Ut 



approaching the limiting value 



lim k = IT = 6. (32) 



This means that the onset of disorder at this particular phase transition is dominated by the creation of 
defects of the fixed width IT ■ In particular, one can check that the width of all the defects in fig. 0(c) 
is about 6 lattice spacings. 

All of the order parameters defined above take a zero value in the non-chaotic phase and increase 
continuously from zero in the chaotic phase as we move away from the transition point. What is more 
interesting, they scale algebraically with the deviation of u from the transition point (fig. 0), all with 
different critical exponents: 

(33) 
(34) 
(35) 
(36) 

One can use these values to compare the transition with similar phase transitions in other conserving 
system, both continuous and discrete, and perhaps determine whether this transition belongs to some 
universality class. 



^max C>C [Uc 


-uY\ 


/3a P 


:i 0.8 ±0.03; 


h OC (Uc - 


-uf'-, 


/3h^ 
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Pt OC (Uc - 


-uf^, 


/3p^ 


s 1.0 ±0.01; 


Cch OC {Uc 


-uf% 
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!i 0.5 ±0.01 



C. Critical exponents 

A closer look at the details of the phase transition at u = Uc reveals that equation (BSh is to be 
expected and that the transition should necessarily be continuous, as a result of the conservation law 
and a particular feature of the local map f{x). 

As was mentioned in the previous section, the change in the growth rate of the k — n mode is 
responsible for the transition. One can determine from eq. (^ that for u close to Uc, equating the 
Lyapunov exponent with the growth rate, the linear stability analysis of the uniform state gives 

A^ « -8b{u - Uc) (37) 

which changes sign as the system moves across the transition point, from the ordered state LI to 
disordered state T2. At the transition point the growth rate obviously vanishes, making the zig-zag 
state neither stable nor unstable in the linear sense. In fact as mentioned in section |l| the zig-zag state 
given by eq. (o) (type-II 2-cycle) with arbitrary amplitude A can exist a.t u = Uc (and only at w = Uc) 
and is stationary, meaning that the amplitude A neither grows nor decays. This is in contrast with the 
result for phase L2 (type-I 2-cycle) , where the amplitude of the stable state is defined by eq. (^ . 

This fact results in some interesting consequences for the system dynamics in the disordered phase 
close to the transition point. Most of the lattice develops a zig-zag pattern (fig. ^ similar to the one we 
just discussed 

<=u, + (-l)^+"A„ (38) 

where now the amplitude Ai is not a constant, but a slowly varying function of the lattice site. The whole 
lattice cannot be in such a state for u < Uc because of the conservation law. In order to compensate 
for the difference, several similar turbulent defects separating laminar domains form, with a fixed width 
I ~ If (see previous section) and the local density of the map variable 

, io+i — 1 

7 Y^ ul' = Uc-5u (39) 

i=io 

which is lower than the critical value Uc by 5u. 

Here 5u, should not strongly depend on u, because u is not a local parameter prescribing the dynamics. 
On the contrary, the local density in the turbulent defect is only determined by the structure of the 
interface separating two laminar domains that have their local densities fixed at Uc, independent of u. 
The structure of the interface, in turn, depends primarily on the width of the turbulent domain, which 
is seen to depend very weakly on u. Numerical results support this conclusion. 

This results in the value of the conserved quantity being "adjusted" to comply with the conservation 
law to give on average 

PiUc + Ptiuc - 5u) = u. (40) 

Now we can easily extract the dependence of pt on u. Since pt + pi — 1 

--'^^ (41, 

This derivation confirms the value of the critical exponent /3 = 1. Thus the conservation law is 
ultimately responsible for the way this particular phase transition occurs and for its type. The ordered 
state turns into a disordered one by developing a set of very similar turbulent domains (defects), which 
have a fixed length (the deviation in eq. (32|) from the It = 6 is due primarily to "defect interaction" 
effects), but whose number increases as the system moves further and further away from the transition 
point, so as to compensate for the change in the density of the map variable in the laminar regions with 
respect to the average value given by the conservation law. We may therefore suggest that defects play a 
more important role in order-disorder transitions in conserving systems than in non-conserving systems 

The value of the critical exponent for the maximal Lyapunov exponent f3\ « 0.8 numerically obtained 
in eq. ( p3| ) is different from the one predicted by eq. (p7|). However the latter value is calculated for a 
reference trajectory corresponding to the uniform configuration and since for u < u^ the system is in the 
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chaotic state, the vahdity of eq. ( p7| ) is far from being obvious, no matter how small the distance to the 
transition point is. We should also consider finite size effects. It is natural to expect that numerical and 
theoretical values agree if there is just one positive Lyapunov exponent and therefore no mode mixing. 
From eq. (|g) it follows that the next mode to become growing is the mode with fc = tt — j^tt and this 
happens when 

u = u,-Au, Au^j^. (42) 

So we might expect the crossover behavior for the exponent /3>: 

/3a ~ 0.8, uc-u> Am, (43) 

/3a = 1.0, Uc-u< Au. (44) 



In practice the value of Au was usually so small (Au « 3 • 10^'^ for L = 128), that eq. (O) was satisfied 
for all deviations of u from the critical value used in our numerical calculations (at most 2 • 10^^). 

And finally we would like to mention that the value of the critical exponent (3h can be evaluated 
from the limiting form of the Lyapunov spectrum. Close to the transition point the positive Lyapunov 
exponents are well approximated by the following expression: 

A,„ = Amaa; -cito'^, 7 = 0.7 ±0.1 (45) 



Together with the equations ( p4D and (33) this implies that the critical exponent corresponding to the 



Kolmogorov-Sinai entropy density is given by 

Ph = /3a (1 + -), (46) 

7 

that yields the value f^h = 1-94 ± 0.20 consistent with eq. (pi]). 

We might expect these exponents to define a universality class for the onset of spatiotemporal chaos. 
The arguments leading to the predictions for the values suggest that the class may depend both on the 
existence of a conservation law and on special symmetry properties of the map function. The restrictions 
on /(x) can be obtained in the following way. 

One starts with the relation between the amplitude and the local density of the type-II 2-cycle for an 
arbitrary f{x): 

E ;7i^""V^"^H = ^ (47) 



n 

n — odd 



Close to the transition point u = Uc it can be rewritten as 

= -r(u - Uc)A + d^A^ + dsA^ + ■■■ (48) 

where r = f"{uc) and coefficients dn are defined as 

dn^-j'^^'Huc), n>3. (49) 

Depending on the values of parameters r and da, we can have either a subcritical or supercritical 
bifurcation at the transition point. The special case studied in this paper corresponds to 

dg - da = • • ■ = (50) 

and is intermediate between the two types of bifurcation. 

These special properties follow immediately for continua of a and b from the parabolic nature of the 
map we have used. More generally we can consider maps of the form 

f{x) — ax + bs{x), (51) 

where s{x) is an arbitrary function symmetric about its maximum. By rescaling and shifting the origin 
of x and choice of normalization of s we can set s(0) = 0, the maximum of s to occur at s = 1/2 and 
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then s(l/2) = 1, leaving the two parameters a and b as well as the conserved quantity u to define the 
system. For this general family of maps the degenerate bifurcation to the period-2 state occurs only for 
a = 1/2 and for u at the maximum of s i.e. u — 1/2. Thus the universality class is codimension-2 — 
two parameters must be tuned to arrive at this type of transition. 



For other values of a and b in (51) the bifurcation to the period-2 state will be either supercritical 
or subcritical. For the supercritical case a stable laminar 2-cycle state develops, which may be the first 
step in a subharmonic cascade. For the subcritical case attractors develop far away in phase space, and 
a full non-linear analysis is needed to determine the type of behavior. 



D. Hysteretic transitions 

Now we turn our attention to the phase transitions that we expect to occur at the boundaries L2-T1, 
L2-T2 and Ll-Tl. As we are going to see later, all three are very similar, so we will concentrate on the 
transition at L2-T1 below. 

There is a considerable difference between the transition at L1-T2 and the transition at L2-T1: in the 
former case the asymptotic state in both LI and T2 is unique, while in the latter case the asymptotic 
state in the ordered phase L2 can be either ordered, or spatiotemporally chaotic. As a result we should 
specify between which states the transition occurs. 

It will be convenient to introduce an additional parameter ujw) characterizing the volume of the of the 



basin of attraction. For example, as we know from section IV in the sub-phase L2i {up < u < Un) the 
attractor is unique and therefore the basin of attraction is the whole configuration space with volume 
vi{u) ^ l^-^ ^ I. 

In the sub-phase L2p the non-chaotic attractor coexists with the chaotic one, so we have < vi{u) < 1 
for Ua < u < Up. Since there are no other attractors in this sub-phase, the volume of the basin of 
attraction of the chaotic attractor is given by vt{u) = 1 — vi{u). Numerical data suggest that most 
of the initial conditions in L2p result in spatiotemporally intermittent chaotic asymptotic state and 
therefore typically Vt{u) ^ vi[u). Moreover viiu) ^ as u -^ m+ while Vt{u) ^ as u ^ u~. 

At the point u = Ua the 2-cycle state loses its stability through a subcritical bifurcation. Outside 
the phase L2 the 2-cycle state cannot exist and as the system crosses the phase boundary all order 
parameters jump from zero in the phase L2 to some nonzero values in the phase Tl (see fig. 0). As a 
result one observes a discontinuous transition in the non-chaotic state: the laminar state abruptly turns 
into the chaotic one. 

However, there is apparently no phase transition in the chaotic state at this point: the dynamics of 
the system on the Tl side of the boundary is very similar to the dynamics of the system evolving on 
the chaotic attractor on the L2 side. All order parameters change continuously as the system crosses 
the boundary from L2 to Tl. Since vi(u) — > as m — > uj, changing the direction does not modify this 
conclusion: the chaotic attractor of the phase Tl smoothly transforms into the chaotic attractor of the 
phase L2. 

Ax u — Up the situation is reversed. Obviously there could be no phase transition in the non-chaotic 
state. On the other hand, the chaotic attractor does not exist in i2;, therefore there should be some 
kind of phase transition in the chaotic state: the order parameters are zero for Up < u < Un, but take 
on nonzero values in the chaotic state for Ua < u < Up. Thus u = Up corresponds to the onset of STL 

The coexisting attractors form a hysteresis loop in the sub-phase L2p (fig. 0). If we start a,t u > Up 
and gradually decrease parameter u the system will remain in the non-chaotic state while u > Ua and 
then jump to the chaotic one at u = Ua- Conversely starting at m < u^ and gradually increasing 
parameter u makes the system remain in the chaotic state while u < Up. At u = Up the chaotic state 
becomes non-chaotic thus closing the loop. 

There is a numerical complication here: it is not possible to establish the exact value of the critical 
parameter Up for a system of finite size. For u -^ Up the lifetime of the chaotic transients becomes 
very long and one cannot reliably determine the type of the asymptotic state. On the other hand, 
Vt (u) — > as M — » Up , which means that it becomes increasingly hard to find initial conditions resulting 
in a persistent chaotic state as u gets close to Up, especially for small systems. Although vt{u) grows 
rapidly with the system size, the smaller \u — Up\ is the larger the system should be in order to obtain 
the persistent chaotic state. Therefore one cannot reliably determine the values of the order parameters 
in the chaotic state close to u = Ur,- 
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As a result, it is even hard to determine reliably whether the transition at u = Up in the chaotic state is 
actually continuous, although the correlation length seems to diverge approaching the transition point. 
The order parameters do not provide a clear picture either. Both h and pt seem to increase continuously 
from zero, but Xmax and Cch jump discontinuously as the system moves across the transition point to 
the chaotic sub-phase L2p. All order parameters gradually increase as the system moves away from the 
transition point toward the chaotic phase Tl. The number of positive Lyapunov exponents also grows 
showing the increase in the number of chaotic modes. These results are consistent with |gj|. But so far 
we do not have reliable data that will allow us to determine whether this transition in fact belongs to 
the universality class of directed percolation or not. 

There is a similarity between the transition at u = Up and the continuous phase transition at u — Uc'- 
the laminar state becomes chaotic through the appearance of a set of turbulent domains that gradually 
spread over the whole system. Nevertheless there is an important difference: the STI state observed 
inside T2 can never become completely laminar because of the conservation law, while the STI state 
inside L2p may decay into a completely laminar state with time. 

The above discussion applies completely to the chaotic sub-phase L2„ which corresponds to u„ < u < 
Ub- One just has to replace Tl with T2, Up with u„ and Ua with Uf,. In particular, vi{u) — > as u — > Uj^" . 
One particular feature of the chaotic attractor in this sub-phase is worth mentioning: for 0.1 ^ w ^ 0.15 
the measure of the turbulent set, pt grows linearly with u. 

The phase transition Ll-Tl at u = u^ is also very similar to the transition L2-T1. One just has to 
replace L2p with L\f, Up with Uf, and Ua with Ud in the above discussion. 

Here, unlike the phase L2, the laminar state is the uniform state, not the 2-cycle. The uniform state 
is unstable for m > u^;, so in order for the laminar regions to be stable the local density inside them 
should be in the range Uc < u < Ud characteristic of the ordered phase LI. This is in fact the case: the 
value of the local density in the laminar regions is ui « 0.53. 

Inside L2/ the number of positive Lyapunov exponents grows linearly with pt , while the latter grows 
linearly with u. This is the kind of behavior one expects in a system where all chaos is localized in 
turbulent domains. 

All three cases of hysteretic transitions observed in our model are very similar and represent many 
of the characteristic features of the specific route to chaos. In particular, STI appears whenever the 
uniform absorbing state experiences a subcritical bifurcation at the transition point (e.g. u = Ua)', the 
hysteresis loop forms as a result of the bistability; onset of STI (e.g at u = Up) can be either continuous 
or discontinuous, depending on the types of defects supported by the laminar (absorbing) state; the 
lifetime of the "laminar" transients diverges at the onset; the confinement effects enhance the stability 
of the absorbing state in relatively small systems. 



VI. LYAPUNOV DIMENSION 

In this section we will briefly comment on the applicability of the Lyapunov dimension as a parameter 
characterizing the dynamics in the system under consideration. This has been used by a number of 
authors (e.g. |^,|l^) to assess the strength of chaos in nonlinear dynamical systems. 

The Lyapunov dimension is defined as 

DL = n+ '^ , (52) 

l^n + 1^71+1 

where z/„ = X]"=i ^i ^^'^ " is such that ;/„ > 0, Vn+i < 0. It was suggested in ref. M that D^ is not 
defined for the system considered, because for the values of the control parameters used all z^„'s were 
positive. 

This conclusion appears incorrect. First of all, this only happens for some restricted set of control 
parameters. For other values of parameters the Lyapunov dimension is perfectly well defined. Second, 
there is no problem with all I'n's being positive: it only means that D^ = L ~ 1, i.e. the dimension of 
the attractor coincides with the dimension of the configuration space (L variables with one constraint) 
meaning that the attractor fills the configuration space. As one can see from the fig. M this only 
happens away from the boundaries inside the chaotic phases Tl and T2, where the system is in the 
strongly chaotic regime with pt « 1. 

The numerical results show that the Lyapunov dimension thus defined is a good (albeit costly to 
calculate) measure of the strength of chaos in the system under consideration. It is a continuous 
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function of control parameters and can in principle be used as an alternative order parameter, though 
it proved to be hard to calculate with the necessary precision close to the phase transition points. 



VII. SINGULARITIES IN THE LYAPUNOV SPECTRUM 
A. Introduction 

It has been suggested 0] that all coupled maps with a conservation law should have a singularity in 
the spectrum of Lyapunov exponents at the value A = 0. The origin of the divergence of the number of 
Lyapunov modes with negative exponents close to zero is generally explained by the following arguments. 

A spatially uniform autonomous system with a locally conserved density is considered. The suffi- 
ciently coarse grained asymptotic state of such a system is supposed to be uniform, i.e the density in 
the asymptotic state should not depend on the spatial coordinates. Now a long wavelength density per- 
turbation is imposed on such an asymptotic state. Because of the conservation law the only mechanism 
by which such a perturbation can decay is diffusion from the regions with high average density to the 
regions with low average density. This process can be described by an effective diffusion equation, which 
can be considered as a long wavelength approximation of the original evolution equation. As a result 
the decay rate of a long wavelength perturbations should be given by a quadratic function of the wave 
vector k. 

In our model this functional dependency can be motivated by analogy with the analytic expression 
obtained for a stable non-chaotic phase LI. Since there is an exact 1-to-l correspondence between Fourier 
mode number and the decay rate in LI, we can formally rewrite eq. (O), re-expressing the number of 
the exponent m through the wave vector k: 

\{k) = \n\l - 4:{a + b - 2bu) sm^{-)\. (53) 

Here for fc — > we have A(fc) oc fc^ and the density of Lyapunov exponents diverges 

n(A) oc lAr^/^ ^ cx), A^O". (54) 

One can attempt to apply a similar numerical analysis in the chaotic phase. We start with writing 
the evolution equation (W) in Fourier space: 



u 



"+!—„," Af^ I k Oh.,\.A„'il \.,n I AU,.\„2/^\ Y^ 



ut-Aia + b-2bu)sm'i-)ul+Absm'i-) ^ ul,ul_^,. (55) 



Consequently the Jacobian Jj"^ = du^^^ / du^ of the evolution transformation takes on the following 
form in the Fourier space: 

J,"fc, = (1 - 4(a + 6 - 2bu) sin2(^)) 4fc' + 86<_fe, sin2(^)(l - 4fc')- (56) 

Fourier modes are not eigenvectors of this matrix (unless all uJJ = as in the uniform state) , but since 
the off-diagonal elements are of order O(fc^) Fourier modes can serve as a good first order approximation 
to the exact eigenvectors for k sufficiently small. This argument suggests that wavevector k should be a 
good label for the slow modes (and small A) of the system in chaotic as well as spatially ordered states. 

B. Structure of the Lyapunov Vectors 

Numerically the Lyapunov spectrum can be calculated using the QR decomposition of the product of 
the Jacobians J/J^^ 

J"Q"-i = Q-'^R^^ (57) 

where Q is orthogonal and R - upper-diagonal to yield: 
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J"... jO = Q"i?"...i?o = Q"^". (58) 

The columns of Q" give Lyapunov vectors and the diagonal elements of i?" - the corresponding Lyapunov 
exponents on the n-th step: 

A„ = i HR^J =< IniR'^J >n . (59) 

Figs. pX|-p^ show several typical time averaged power spectra of the instantaneous Lyapunov vectors 
along with the corresponding Lyapunov spectra. The power spectra are represented in the form of the 
density plots showing the relative contribution R„i{k) from the Fourier mode with number kL/iii to the 
?7i-th Lyapunov vector, while the Lyapunov spectra show the correspondence between Lyapunov vectors 
and exponents. The power spectra are normalized 

^F,„(fc) = l, (60) 

k 

so that Rm{k) can alternatively be interpreted as the probability distribution functions. 

We would like to mention that due to the ergodicity of the dynamics the power spectra (as well as 
Lyapunov spectra) should be reproducible, i.e. different initial configuration with the given set of control 
parameters should produce the same spectrum and this is found to be the case. It is interesting to note 
though, that we obtained a unique form for the spectra even for the values of the control parameters 
producing frozen patterns, where we expect ergodicity to break down. 

We start with the spectra calculated inside the chaotic phases Tl and T2. One can easily notice that a 
singularity appears in the spectrum of Lyapunov exponents (fig. QO(b)) when the dominant contribution 
to the Lyapunov vectors, corresponding to a slow evolution (small A) comes from the long wavelength 
Fourier modes and there is at least an approximate 1-to-l correspondence between the exponent number 
and the dominant Fourier mode number for the relevant vector (fig. n^(a)). It is therefore natural to 
expect the small negative Lyapunov exponents to be determined by the decay rate of the corresponding 
Fourier mode. 

In the long- wavelength limit fc — > 0, the latter could be determined from a hydrodynamic analysis [|5[ 
of the problem. It was argued m that at long wavelengths the effect of all short wavelength corrections 
to the equation of motion ( p5[ ) can be combined into a stochastic noise term and a renormalized diffusion 
constant D, producing an effective (discrete) Langevin equation: 



u 



k 



- -Dk'ul + kWk (61) 



with ^-correlated noise term 



<r?rr?? >=B5u'5nn-. (62) 



According to the Central Limit Theorem such noise averages to small values on large lengthscales and 
therefore, for small fc, the decay rate of the Fourier mode u^ is equal to —Dk^. 

Because of the conservation law one of the exponents in our model is always equal to zero. Let it be 
Xmo- It obviously corresponds to the Lyapunov vector with fc = 0. Typically, for small A, the relation 
between the number m of the Lyapunov vector and the dominant long-wavelength Fourier mode k can 
be represented by the following simple form: 

km = a—{m -mQ). (63) 

Parameters mo and a ~ 0.5 are determined from the Lyapunov spectrum and power spectrum, respec- 
tively (mo — 0,a — 0.5 inside LI). 

The Lyapunov exponent X,n is in fact calculated as a time averaged decay rate of the corresponding 
instantaneous Lyapunov vector. Since the time averaged power spectrum Pm(fc) gives the averaged 
relative contribution of the Fourier mode k to the m-th Lyapunov vector, we might estimate the value 
of the TO-th exponent as 

X,n = -DY,Pm{k)k\ (64) 
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As one can deduce from fig. |lO|(a), for u = 0.7, Pm{k) is sharply peaked at km, so approximating 
Pm{k) = 6k^k we readily obtain that the form of the Lyapunov spectrum for small negative values of 
the exponent should be given by 

Xm w -Dkfni m = m-o, mo + 1, . . . (65) 

One can easily check that the numerically calculated exponents are in fact given quite precisely by this 
expression with D « 0.32 for < A; ^ 7r/4. So we recover eq. (|54|), but now for the chaotic state. 

Fig. O provides us with another typical example of the spectra corresponding to the strongly chaotic 
dynamics, now calculated for u — 0.3. One can notice that the singularity in Lyapunov spectrum in fig. 
0(b) is still present, although the quadratic fit provided by ( pq ) is quite poor compared to the one for 
w = 0.7. 

The power spectrum (fig. 0(a)) shows that, similarly to the previous case, the dominant contribution 
to the Lyapunov vectors corresponding to slow evolution comes from the long-wavelength modes, but 
now the distribution Pm{k) is much broader. This means that we can no longer approximate Pm{k) by 
5k,-nk and there could be considerable corrections to eq. (pq), which does not however change the general 
conclusion about the presence of the singularity in the spectrum of Lyapunov exponents at A « 0. 

These arguments work well in the area of strong chaos, where the large scale dynamics is determined 
by the long wavelength modes, i.e. modes with k <^ 1/^, where ^ is the correlation length in the system. 
When we approach a continuous phase transition, though, the correlation length grows and becomes 
comparable to the system size L. When this happens, the nonlinear terms in eq. (|55| ) become relevant 
on every scale [nl and their effect can no longer be emulated by an effective noise term. This results 
in a strong coupling between different long wavelength modes (i.e. modes with k ^ 1/i) and therefore 
the approximate 1-to-l correspondence between the mode number and its growth rate is completely 
destroyed. This means that there is no reason to expect the divergence of the number of slowly evolving 
modes anymore. 

Fig. |lj(a) suggests that there is a strong mode mixing in the system and it is not possible to extract 
the dominant contribution to the Lyapunov vectors. In fact, we really have a continuous phase transition 
at u — Uc, very close to u — 0.455, and it is responsible for the disappearance of a singularity in the 
spectrum in fig. |l2|(b). One might still hope to find a trace of a singularity at this particular value of u 
by going to much larger system size (higher resolution). 

C. Positive side singularity 

Some models with a conservation law [Bj are found to have a density of exponents diverging on the 
positive A side 

n(A)^oo, A^0+, (66) 

although this feature is not considered very common. This divergence is present in our model too (see 
fig. |l0|(b) for example), but is somewhat weaker than the divergence at negative A: for some values of 
the parameters it can be barely seen even for large lattices (typically L > 10^ is necessary). For A — > 0^ 
the spectrum can usually be fitted quite precisely by a quadratic function 

Am ~ D{m — toq)^, m = mo, mp — 1, . . . (67) 

which means that the singularity of n{X) is inverse-square-root on both positive and negative sides. 

Numerically calculated spectra show no sign of smoothing out of the singularity (on either side) with 
increasing resolution (increasing size of the system) up to L = 512 (the results for L = 128 and L = 512 
are presented in fig. Il^(a)) and also suggest that the fit ( |67| ) could be good for positive A as large as 

Figs. |l^(a) and |ll](a) suggest that it is possible to extract the dominant Fourier modes corresponding 
to small positive A. Again their wavevectors scale roughly linearly with toq — m and therefore equation 
( pSi) should be replaced by 

km = a—{mo-m) (68) 

for m < mg. Equation (^) though should be abandoned in favor of a more precise one, preferably 
derived directly from the evolution equation (pq). We intend to explore this question in more detail 
later. 
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D. Effective diffusion constant 

We should also mention that the values of the effective diffusion constant D numerically obtained 
from eq. ( pq ) do not coincide with the ones obtained through the dynamic structure function S{k,t). 
According to (|6^) it is defined (see [Q) as 

S{k,t) = Bk^e-''''"'. (69) 

It will be convenient to distinguish these using the notation Dis for the former and Dgf for the latter. In 
the regions of parameter space where the effective Langevin equation is applicable (and hence eq. ( [69| ) 
is valid) these approaches can give substantially different results. This is yet another indication that 
equation (p4) can only provide a very crude estimate and should be amended considerably to obtain 
adequate results. 

It is instructive to compare how Dig and D^f change if parameter b is varied while both a and u are 
fixed (see Fig. |3|). For u > 1.3 Di^ > Dsf and they both grow with increasing b. Numerical data 
available so far suggests that Dis/Dgf -^ 1 as 6 — > oo (for strongly chaotic systems). 

For w ^ 1.3 Dsf drops almost to zero indicating that there is almost no diffusion in the system. Indeed 
we know that this is the region of locked chaotic dynamics (Tli). So this result is not surprising: the 
formation of a locked structure prevents diffusion (on any scale larger than some typical scale determined 
by that structure). We cannot probe smaller scales using the dynamic structure function because the 
effective Langevin equation l pi\ ) is not valid for k large enough, but supposedly diffusion survives there. 
It is interesting to note however that locking has apparently no effect on the Lyapunov spectrum: the 
change in Dig is very gradual across Tl. 



VIII. ROLE OF THE CONSERVATION LAW 

Now that we have studied the dynamics of the CML with the conservation law quite thoroughly and 
compared the characteristic phenomena with those observed in CML's without any conservation laws 
we would like to discuss whether the distinguishing features are really explained by the presence of the 
conservation law. 

It is often very hard to distinguish the effect of the conservation law on the dynamics of a system from 
the effects introduced by other aspects of the evolution equation. Here we would like to explore some 
consequences of violating the conservation law, trying to retain the structure of the original equation 

First of all, we are looking for internal homogeneous perturbations that would violate the conservation 
law of eq. (M, but would preserve the structure of the equation. This is easily furnished by the following 
modification of the original evolution equation: 



«"+' = u^' 



eg{u-) + {f{uU) - 2/«) + /«+i)), (70) 



where the local map function ,f{x) is the same as above and a perturbation g{x) is introduced. One can 
easily notice that in the simplest case of the uniform state it reduces to the equation determining the 
evolution of the average density 

u"+i=u" + e5(u"). (71) 

This equation, though, does not provide us with any reliable information concerning the dynamics of 
the average density in the case of a nonuniform state. 

Second, in order to be able to compare two systems, one with and one without the conservation law, 
we should ensure that the latter is violated only "mildly" . In other words, we would like the perturbed 
system to have a phase diagram that could be compared to that of the original system. In particular 
we would like to preserve the dimensionality of the parameter space. Since u is no longer conserved 
it is not a parameter of the dynamics. Instead we introduce another parameter uq that will enter the 
evolution equation through the perturbation function g{x). 

By "mild" conservation violation we also mean that the change of the originally conserved average 
density u during any single time-step is sufficiently small: 

|m"+1-u"|<1. (72) 

17 



We would also like the fluctuation of the average density u to be bounded, such that 

uq ~ Su < m" < uq + 6u (73) 

for some finite 6u at any time step n. Then we would be able to compare the dynamics of the perturbed 
system with that of the original system with conserved quantity u « mq- 

Our primary interest in this section is the relation of the conservation law to the existence of a 
singularity in the spectrum of Lyapunov exponents. It is thus reasonable to compare the perturbed 
system with the original one, described by the equation (n|), at the value of the conserved quantity 
u — 0.8, where the conserving system displays a strongly chaotic dynamics and has a Lyapunov spectrum 
with a pronounced singularity at A = 0. 

We will start with the following choice of the (nonlinear) perturbation function: 

g{x) - («o - 2;)3 (74) 

and choose e = 0.2. The dynamics of the perturbed system seems qualitatively very similar to the 
dynamics of the unperturbed system, but now the average density u is not conserved and fluctuates 
about Mn with a standard deviation of order few percent. 

Fig. |lj of the Lyapunov exponent spectra focuses on the parts corresponding to slow evolution (small 
A). Comparing the spectrum of the modified system with that of the original system for u = uq = 0.8, 
one can easily notice that the singularity is clearly present in the conserving case, while the spectrum 
of the perturbed system appears to be similar to the spectrum of the original system, but slightly tilted 
and shifted downwards. 

In order to better understand the origin of such a metamorphosis it is advantageous to use another 
type of perturbation which is a lot easier to interpret and study analytically: 

g{x) =uo -X. (75) 

Since g{x) is linear, equation ( |7l| ) now describes the evolution of the average density u for arbitrary 
initial state. For e > (0.001 < e < 0.1 was used) the asymptotic state is a configuration with u = uq. 
The effective Langevin equation corresponding to ([75|) should read 



dtu{x, t) = -eu + Ddlu + dlri{x, t) (76) 

On large lengthscales noise averages out and we obtain the following dependence of the growth rate 
on the wave vector: 

X{k) = -e-Dk^, (77) 

i.e. we should expect (for small fc's) the Lyapunov spectrum of the perturbed system to be shifted 
downwards by i5A = — e with respect to the spectrum of the conserving system, while retaining the same 
type of singularity. The total decay rate is then determined by a linear combination (at least for small 
enough coupling) of the diffusion with local dissipation. 

Comparin g th e spectra of Lyapunov exponents of the modified system (fig. 03) and the original 
system (fig. |lj(a)) for u = uq = 0.8 we see that the numerically obtained spectra of the perturbed 
system follow the prediction of the Langevin equation (UW quite precisely for small (e ~ 0.01) as well 
as for relatively strong (e = 0.1) perturbations. 

We expect the negative shift of the "slow" part of the spectrum to be attributed to the dissipative 
nature of the perturbations used above. In fact we could have a positive shift or no shift at all. In 
order to see this we pick the nonlinear perturbation function g{x) with the first derivative which is not 
negat i ve- definit e : 

{{uo - xy + t) 

with e > 0. We used 7^ = 0.001 and 0.01 < e < 5. 

Equation ( |7l| ) does not hold anymore, but the numerical data suggests that the dynamics of the system 
is ergodic and the averaged density u in the asymptotic state fluctuates about uq with fluctuations being 
again of order few percent. 
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Fig. |l^ shows the Lyapunov spectra of the perturbed system with uq = 0.8. This figure suggests that 
the spectrum is indeed shifted in the direction of positive rather than negative values of A. A small 
perturbation (e = 0.1) results in a small distortion of the original spectrum: the slope of the spectrum at 
the value A = on both positive and negative sides becomes nonzero, i.e. the singularity in the density 
of Lyapunov exponents disappears. The slope increases with e and for sufficiently strong perturbation 
(e > 2.0) all traces of the singularity vanish. 

These examples suggest that there are, in fact, two different aspects of the singularity in the Lyapunov 
spectrum of a system with the conservation law. The first one is the presence of a singularity at some 
value A = Aq. The numerical data obtained suggests that the singularity survives in the special case 
of the linear perturbation function g{x). A nonlinear perturbation results in the disappearance of the 
singularity. The Jacobian of a perturbed system can be written as 

Jkk' = Jkk' + E 7-^ n E <^(E ^' - ^ + fc'), (79) 

where J^^,, is the Jacobian dS^) of the conserving system. If g{x) is nonlinear the off-diagonal elements 
of J^j,, become of order 0(1) instead of 0{k'^). As a result Fourier modes are no longer good as an 
approximation to the exact eigenvectors even for small k. The wavevector k can no longer label the 
slow modes of the system and therefore there is no reason to expect the singularity in the density of 
Lyapunov exponents to remain. It becomes "smoothed out" by the perturbation. 

The second aspect is the actual value of Aq in case a singularity is present. We already saw that 
imposing linear perturbation (^5|) made Aq become negative (cf. (|77|)). This is a consequence of the 
shift of the "slow" part of the spectrum as a whole in response to some local effects, e.g. dissipation. 

We may therefore suggest, that since all known CML's with an additive conserved quantity possess a 
spectrum of Lyapunov exponents distinguished by the presence of the singularity at A = 0, a conservation 
law is a sufficient condition for the existence of such a singularity. It seems to be a necessary condition 
as well, at least in the class of models (|^) studied in this paper. The singularity at A 7^ indicates that 
there is a mechanism of (local) dissipation concurrent with diffusion and if this is eliminated the system 
becomes strictly conserving. This refinement might be helpful when looking for hidden conservation 
laws using Lyapunov spectra [BJ, numeric or experimental. 

All three types of perturbation studied above are seen to have an effect on the phase diagram. For e = 
phase diagrams in 3-dimensional parameter spaces {a,b,u) and (a,5, uq) obviously coincide. Numerical 
data for all types of perturbation studied suggest that gradual increase of the parameter e makes the 
phase diagram of the perturbed system change continuously. The boundaries of the phases shift, making 
some of the phases shrink or completely disappear (e.g. sufficiently strong perturbation of any type 
obliterates the phase L2). Other phases may expand as the stability of their basic state is enhanced 
(as is the case for the phase LI and the uniform state in the presence of a "dissipative" perturbation). 
We can say that, although the phase diagram is sensitive to the violations of the conservation law, it is 
robust with respect to sufficiently weak violations. 

Since the changes in the phase diagram provoked by perturbations are found to be continuous, it is 
reasonable to expect that the phase transitions which were 1-order in the conserving system will remain 
1-order after a perturbation is imposed. This leaves us an opportunity to observe how the characteristic 
exponents change with increasing perturbation. In our model the continuous phase transition at the 
boundary T2-L1 seems to be a promising point of investigation. Numerical data we have at the moment 
does not allow us to answer an important question of whether the universality classes change if the 
conservation law is violated, but we plan to return and investigate this later in more detail. 



IX. CONCLUSIONS 

We have systematically investigated the properties of a coupled map lattice with dynamics constructed 
to satisfy a conservation law and to show spatiotemporal chaos. 

The conserved quantity provides an additional control parameter: as its value is changed, a rich 
phase diagram with a number of phase transitions between ordered and disordered states is found. Both 
continuous and discontinuous transitions occur, as in coupled map lattices without a conservation law. 
The basic structure of the phase diagram is given by the linear stability boundaries of the ordered 
phases, although near the discontinuous transitions bistability may occur. Increasing the nonlinearity, 
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determined by the parameter b (see fig. 1(b)), renders the spatially uniform ordered phase (phase LI) 
unstable. For < u < 1/2 the linear instability of the uniform state occurs via spatial period doubling 
(the zone boundary mode goes unstable). For a < u < 1/2 the transition is immediately to a chaotic 
state, which takes the form of an increasing number of "turbulent" regions of roughly fixed size (defects) 
moving through the "laminar" background. This transition is continuous with a diverging correlation 
length and other scaling approaching the ordered state. For < u < a the transition to chaos is 
through a series of two subharmonic bifurcations, passing first to an intermediate 2-cycle state (the L2 
phase). The onset of chaos from L2 is hysteretic. A complete subharmonic cascade is not observed. For 
1/2 < M < 1 the instability of the LI phase occurs through modes of all wave vectors going unstable 
together, and the appearance of chaos is also hysteretic with frozen chaotic domains developing. 

As in thermodynamic systems, the phase transitions can be conveniently described by the use of order 
parameters, although since we are concerned with the growth of disorder the choice of the appropriate 
order parameter here is by no means obvious. In the case of the continuous transition we find that 
a number of proposed order parameters scale with distance to the transition point, thus allowing the 
evaluation of critical exponents, which may help to pin down whether universality classes for the onset of 
spatiotemporal chaos exist. An interesting question is whether the conservation law, which clearly affects 
the dynamic correlations, also changes the exponents and the universality class (if such a classification 
exists) of the phase transitions. The transition from LI to r2 in the conserving model and in perturbed 
versions where the conservation law is weakly violated should provide a good arena for investigating 
this. 

A symbolic description of the dynamics, reducing the complex states to regions of laminar (ordered) 
and turbulent (fluctuating) regions is useful in describing the chaotic states near the transitions. In 
particular we see that the onset of chaos always happens in the form of turbulent regions gradually 
spreading over the laminar background. Whenever the chaos appears to grow continuously, the turbulent 
fluctuations appear in the form of turbulent defects, usually propagating across the system, with size 
It of order a few lattice spacings. Similar results were obtained by Kaneko |l2|| for the non-conserving 
CML, although in our case there seems to be a stronger tendency for the defects to propagate with 
constant velocity. 

We have also studied the Lyapunov eigenvalues and eigenvectors of the chaotic states. A conspicuous 
feature of conserving models is the singularity in the density of Lyapunov exponents around A = 0. 
There is growing evidence that the singularity is associated with the existence of Lyapunov eigenvectors 
that are labelled by the wave vector k for small k — the singularity is then the usual Van Hove singularity 
coming from mode counting and the assumption of a smooth spectrum at small fc, i.e. A = Aq — D^ffk^. 
This leads immediately to a |A — XqI^-^^^ singularity in the density of exponents. 

Some evidence for this idea comes from the work of Bohr et al. who studied some models where 
the Lyapunov spectrum can be exactly calculated, and in which the eigenvectors are trivially Fourier 
modes, since the Jacobian is independent of the dynamic variables. In these models the singularity can 
occur at nonzero Aq, and the singularity may be to larger or smaller A. 

For the conserving (and nontrivial) maps we have studied, the labelling of the Lyapunov eigenvec- 
tors by the wave vector and the expression for the Lyapunov spectrum, with Aq = and D positive, 
is suggested by a hydrodynamic analysis, which associates the Lyapunov eigenvectors with the diffu- 
sively decaying modes given by a long wavelength Langevin description. The Fourier power spectrum 
of the Lyapunov eigenvectors provides support for this explanation: for small negative A the power 
spectrum corresponding to the m-th eigenvalue Xm is indeed strongly peaked around a small k with 
k (X {m — mo)2Tr/L with Xjn„ — 0. Note however that the diffusion constant estimated form the Lya- 
punov spectrum may be significantly different from other estimates, and the proportionality constant a 
is greater than 0.5 which would be the value from simply counting the Fourier modes. Indeed, there is 
also a concentration of spectral power towards small k on the positive A side, with a peak wave vector 
which again appears to scale linearly with niQ — m for m sufficiently close to toq. Hence the long wave- 
length Fourier modes contribute appreciably to the Lyapunov vectors corresponding to small positive as 
well as negative exponents, which explains the deviation of a from 0.5. There is no understanding from 
the hydrodynamic approach of the positive eigenvalue long- wavelength modes and associated singularity 
in the density of exponents. 

Adding the linear perturbation (uq) makes Ao nonzero, but k remains a good label. Nonlinear pertur- 
bations that eliminate the conservation law however destroy the singularity, and presumably in this case 
k is no longer a good label. Thus the conservation law appears to be a sufficient condition for Lyapunov 
spectrum singularity, but a complete quantitative understanding of this association remains lacking. 
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FIG. 1. Cross-sections of the phase diagram in a 3-dimensional parameter space (a) u vs. a at b = 1.3; (b) 
M vs. 6 at a = 0.4. The sohd hne corresponds to a continuous phase transition. The dashed and the dotted 
hues denote the phase boundaries where discontinuous phase transitions between chaotic and non-chaotic states 
occur. 
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FIG. 2. 1-dimensional cross-section of the phase diagram. Measure of the set of turbulent domains, pt, is 
plotted schematically (actual numeric results are presented in fig. bl) as a function of parameter u. Solid lines 
correspond to the values for the linearly stable asymptotic states and dashed lines correspond to the unstable 
asymptotic states. Two different linearly stable states coexist in the sub-phases L2p, L2„, LI/. A locked chaotic 
state always forms in T\i and for some initial states in LI/ as well. 



22 



256 




(a) 



256 



32 64 96 128 160 192 224 256 
Lattice site 




32 64 96 128 160 192 224 256 

(C) Lattice site 





224- 


^m 1 1 1 1 




192- 


fli ^\. 




160- 


9 ^^l 


V 




B ^Vflb 


fci 


128- 


^L. ^^B H 


H 




^^^ ^H B 




96- 


K^v Iw 




64- 


H ^v IB 




32- 


n n8 









(b) 



256- 



32 64 96 128 160 192 224 256 

Lattice site 




96 128 160 

Lattice site 



256 



FIG. 3. Symbolic representation of the system dynamics close to the phase transition points; turbulent sites 
are marked black and laminar sites are white, (a) "percolating" state at w = 0.06; (b) "nuclear" state at u = 0.09; 
(c) defect dominated state at m = 0.46; (d) frozen pattern at m = 0.66. Lattice size is 256. 
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FIG. 4. Domain length probability distribution (+ - laminar Pi{l), • - turbulent Pt{l)) in (a) "percolating" 
state at M = 0.06, (b) "nuclear" state at u = 0.09, (c) defect dominated state at u = 0.46, (d) frozen pattern at 
u = 0.66. Lattice size is 256. 
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FIG. 5. Correlation length diverges near the continuous phase transition 
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FIG. 6. Critical behavior of average domain lengths: (a) length of laminar domains diverges algebraically, (b) 
length of turbulent domains converges exponentially 
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FIG. 7. Scaling of order parameters at the continuous phase transition calculated on the lattice with 2048 sites: 
(a) maximal Lyapunov exponent Xmax, (b) Kolmogorov-Sinai entropy density h, (c) measure of the turbulent 
set pt, (d) intensity of chaotic modes Ech 
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FIG. 8. Lyapunov dimension Dl (•) and turbulent set measure pt (+) as functions of the conserved quantity 
u, calculated on the lattices with 128 and 2048 sites respectively. Only the values in the chaotic state are shown 
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0.8- 








a 


0.6- 








1 

03 

> 


•* 


•• 


- 


ft 


0.4- 






_ 


§ 


0.2- 


t 


• 










- 




0- 


1 1 1 


1 1 r 


1 



32 



64 



96 128 160 

Lattice site 



192 224 



256 
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domains with difi'erent amplitudes Ai and A2 are separated by two turbulent defects of the same width. Lattice 
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FIG. 10. Power spectrum (a) and Lyapunov spectrum (b) at it = 0.7, inside the chaotic phase Tl. The 
dominant contribution to the Lyapunov vectors corresponding to small A comes from the long-wavelength Fourier 
modes only (the white dashed line gives the fit provided by eq. (m) with a — 0.6). As a result a pronounced 
singularity appears in the Lyapunov spectrum (exponent density becomes singular at A = 0). Lattice size is 128. 
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FIG. 11. Power spectrum (a) and Lyapunov spectrum (b) at u = 0.3, deep inside the chaotic phase T2. 
The dominant contribution to the Lyapunov vectors corresponding to small A comes from the long-wavelength 
Fourier modes, but due to the contribution from mid-wavelength modes the singularity in the exponent density 
becomes much weaker. Lattice size is 128. 
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FIG. 12. Power spectrum (a) and Lyapunov spectrum (b) at m = 0.455, close to the point where the system 
experiences the continuous phase transition. There is a considerable contribution from the short-wavelength 
Fourier modes to the Lyapunov vectors corresponding to small A resulting in the disappearance of the singularity. 
Lattice size is 128. 
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FIG. 13. Effective diffusion constant as a function of parameter b. a — 0.4 and u = 0.8. The values of Dis 
(+) are determined using the Lyapunov spectrum of the lattice with L — 128, assuming a = 0.5. Dgf (•) is 
obtained from the dynamic structure function calculated for k — O.OItt on the lattice with L = 2048. 
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FIG. 14. "Slow" part ol Lyapunov spectra of the unperturbed system and the system with nonhnear "dissi- 
pative" perturbation, (b) The spectrum of the perturbed system (mq — 0.8, L = 128) with e — 0.2 appears to 
be shifted downwards and tilted with respect to (a) the spectrum of the unperturbed system {u — 0.8). Here 
the data for two lattice sizes, L — 128 (•) and L = 512 (+), is superimposed to show the finite-size effects. The 
density of exponents is seen to be singular on both positive and negative side. The solid and the dashed lines 
give the quadratic fit provided by eqs. (|67[) and (pSl) respectively. 



95 



c< 



(a) 




e< -0.1 



75 85 95 

Lyapunov vector number 



105 




(b) 



75 85 95 105 

Lyapunov vector number 



FIG. 15. "Slow" part of Lyapunov spectra of the system with linear "dissipative" perturbation. For small 
perturbation (a) e = 0.01 as well as for large perturbation (b) e — 0.1 the spectrum gets shifted downwards 
by 5X ft; — e with respect to the spectrum of the unperturbed system. Lattice size is 128 and uo = 0.8. The 
dashed line gives the quadratic fit provided by eq. ([77|) with D = 0.33 (assuming that k is given by eq. ( |63[ ) 
with a — 0.5). 
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FIG. 16. "Slow" part of Lyapunov spectra of the system with "mixing" perturbation. For small perturbation 
(a) e = 0.1 the tangent line (dashed) at the inflection becomes tilted. The singularity disappears. For very 
strong perturbation (b) e — 5.0 the inflection vanishes completely. Lattice size is 128 and uq = 0.8. 
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